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Study of the order of the phase transition in pure U(l) gauge theory with 
Villain action 
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We address the question of the order of the deconfinement phase transition of four dimensional U(l) lattice 
gauge theory. Simulations of the Z-gauge theory dual to the Villain action on toroidal lattices up to lattice sizes 
of 28 4 give results consistent with both, a vanishing and a nonvanishing discontinuity in the thermodynamic limit. 
A decision on the order of the phase transition requires still larger lattice sizes. 
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1. INTRODUCTION 

U(l) lattice gauge theory with Villain action 

exp(-S v ) = n E exp[-^p-27rm P ) 2 ](l) 

P mpE 1 ^ 



is, on the infinite lattice, equivalent 
gauge theory 



i 

w 



to the Z 



(2) 



a gauge theory with an integer valued gauge field 
on the links. In particular the order of the phase 
transition and universal quantities are the same. 
We choose the Z gauge theory for computational 
advantages and use a toroidal lattice with peri- 
odic boundary conditions. We expect[] the Villain 
action to have properties similar to those of the 
extended action 



Se = - 53(/3 cos(6» P ) + 7 cos(2(9 P )) 



(3) 



with 7 k, -0.22. 



The order of the phase transition in this pa- 
rameter region is controversial. Scaling of gauge- 
balls (2| suggests second order. Finite size scal- 
ing of the specific heat is consistent with first 
order J3U3] (torus and surface of 5D hypercube) 



1 This is suggeted by a Fourier transform of S\ 



as well as with second order |5||| (homogenous 
sphere). The volume dependence of the latent 
heat on finite lattices is consistent with first or- 
der as well as with second order. 

For a more detailed discussion of our results we 
refer to [Q. 

2. SIMULATION 

A Metropolis algorithm is used for the up- 
date. The statistics is summarized in Tab. [I]. 
We perform Ferrenberg Swendsen multihistogram 
reweighting Q to the maximum of the specific 
heat. Jackknife errors are calculated for all quan- 
tities. The flip autocorrelation time Tfli p , which 
corresponds to the mean life time of a pure phase 
configuration, is estimated from the time evolu- 
tion of the energy density. Tflj p is a more adequate 
measure for the number of statistically indepen- 
dent measurements than n n t. Tfl; p is in our case 
about a factor 3-5 greater than r- mt . 

3. SPECIFIC HEAT 

Finite size scaling theory for first order tran- 
sitions (see e.g. ||) predicts the maximum of 
the specific heat to diverge oc L D which implies 
v = 1/4. The effective exponent u e ^{L), which 
is calculated from the specific heat at two neigh- 
bouring L, is shown in Fig. [l| We compare with 
v c e obtained in at 7 = —0.2 and find agree- 
ment almost within error bars. 
Clearly there is a tendency for v e s{L) to decrease 
over the L range considered, but a stabilization 
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0.77716 


0.778332(24) 


10.5 


233 






0.77822 




14.1 


228 




10 4 


0.77713 


0.777556(21) 
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0.77706 


0.777183(12) 
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11.5 
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14 4 


0.77699 


0.776976(07) 


9.9 


1697 


9000 


16 4 


0.77689 


0.776882(11) 


7.8 


2852 
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0.77681 


0.776816(07) 


4.42 


4384 
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5156 


19000 
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0.77677 
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8639 
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1.05 
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0.776760(06) 


2.55 
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2.862 


11383 


47000 


24 4 


0.77674 


0.776742(05) 
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16360 
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2.1 


17807 


50000 


26 4 


0.776745 


0.776758(10) 


1.42 


35151 


110000 


28 4 


0.77672 


0.776726(07) 


0.95 


28701 


120000 



Table 1 

Lattice sizes, k values and positions of the spe- 
cific heat maxima. The statistics is given in 10 6 
sweeps. Also included are the integrated autocor- 
relation time and an estimate of the flip duration 
in sweeps. 

on a certain value is not observed. 

4. PROBABILITY DISTRIBUTION 

The probability distribution pi,{e), where e is 
the action divided by the number of plaquettes, 
has for L > 10 two peaks with a minimum in 
between them. It is not well described by the 
sum of two Gaussians. The minimum becomes 
more pronounced with increasing volume, but the 
distance between the peaks decreases. The his- 
tograms are reweighted to the maximum of the 
specific heat. For a first order transition this is 
asymptotically equivalent to equal weight of the 
coexisting phases. Then the positions in e of the 
maxima, the minimum and where the histogram 
is half of its respective maximum value on the 
outer flanks of the distribution function are cal- 
culated. 

A first order transition would obviously require 
the width of pz,(e) to extrapolate to a non-zero 
value. It is consistent with the data that the his- 
togram width shrinks to zero approximately like 
1/L for L — ► oo (Fig. |J). This behaviour would 
imply a second order phase transition in spite of 
the observation of hysteresis and a double peaked 
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Figure 1. The effective exponent u e s(L) for the 
Villain action (|l|) (squares) and the extended ac- 
tion (H) (triangles) at 7 = —0.2. The values for 
the extended action are taken from Fig. 15 in Q|. 

histogram on finite lattices. 

The position especially of the lower maximum 
(confinement phase) varies irregularly. Presum- 
ably the discrete nature of the considered action 
contributes to this effect. 

The L extrapolation of the upper energy 
branch's data (Coulomb phase) is less sensitive 
to the extrapolation ansatz than the lower en- 
ergy branch's data. This is related to the fact 
that the energy fluctuations within the confine- 
ment branch are much larger than in the Coulomb 
phase. In fact they have the magnitude of the 
possible gap itself. 

Following we use 



F(L) = 



1 



2L d ~ l 



log 



pinax,l pmax,2 
L L 



(4) 



as a finite volume estimator for the interface ten- 
sion. F(L) slowly decreases with the lattice size. 
Again a stabilization to a constant value is not 
observed. It may shrink to zero for L — > 00. A 
power law fit F{L) oc L~ z results in z = 0.3(1). 

Also the ratio of the peak heights is not stable 
in the considered L range. 
5. CONCLUSIONS 

We have done a high statistics simulation of the 
Z gauge theory using large lattices. Nevertheless, 
it is not possible to give definite answers to the 
question of the order of the phase transition from 
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Figure 2. The positions of the maxima (squares), 
of the minimum (circles) of the energy density 
probability distribution, and the position on the 
flanks, where the probability distribution has half 
the value of the respective maximum (triangles). 

our simulations. We find that the typical indica- 
tions of a first order transition — hysteresis and 
histograms with double peaks — persist up to 
L = 28. On the other hand, the width of the 
energy density histogram might shrink to zero, 
which would imply second order. The L depen- 
dent exponent v e g{L) decreases steadily but does 
not stabilize in the considered L range. The limit 
limi^oo v $(L) requires larger lattices. The finite 
volume estimator for the interface tension slowly 
decreases with L and might vanish in the thermo- 
dynamic limit. In addition, the flattening of the 
probability distribution in the metastable phase 
region typical for first order phase transitions has 
not been observed in our simulations. Therefore a 
stabilization of the finite volume interface tension 
estimator is not expected. The ratio of the peak 
heights at the maximum of Cy does not stabilize 
in the considered L range. 

The simulations have been done on toroidal 
lattices as opposed to spherical lattices with ho- 
mogenously distributed curvature (||| 7 where full 
consistency with second order finite size scaling 
behaviour has been seen. Also, we want to point 
out that a mass determination on toroidal lat- 
tices away from the transition point — where fi- 
nite size effects are under control — , as has been 
done in is feasible in spite of the complica- 
tions immediately at the transition point. Those 



results favor second order. Our data for the effec- 
tive v e g(L) indicate that if the phase transition 
should actually turn out to be first order, the next 
larger lattice sizes may exhibit a stabilization of 
v e s(L) at the value 1/4. Such lattices could, in 
principle, be simulated on massively parallel ar- 
chitectures. 

The transition might well be of first order, but 
even if that was a given fact larger lattices would 
be necessary for a better understanding of the pe- 
culiar finite size effects. In that case the similarity 
to second order would indicate a critical point not 
far away in the coupling space. In any case it will 
be interesting to pursue the investigation of pure 
U(l) gauge theory on the lattice further. 
6. ACKNOWLEDGEMENTS 

This work is supported in part by funds pro- 
vided by the U.S. Department of Energy (D.O.E.) 
under cooperative research agreement DE-FC02- 
94ER40818. The computations have been per- 
formed on the Cray-T90 of HLRZ Juelich. We 
thank I. Campos, J. Jersak and U.-J. Wiese for 
discussions. 

REFERENCES 

1. M. E. Peskin, Ann. Phys. 113 (1978) 122. 

2. J. Cox, W. Franzki, J. Jersak, C. B. Lang, 
T. Neuhaus and P. W. Stephenson, Nucl. 
Phys. B499 (1997) 371; J. Cox, J. Jersak, 
H. Pfeiffer, T. Neuhaus, P. W. Stephenson 
and A. Seyfried, [hep-lat/9808049| 

3. I. Campos, A. Cruz and A. Tarancon, Phys. 
Lett. 424B (1998) 328. 

4. I. Campos, A. Cruz and A. Tarancon, Nucl. 
Phys. B528 (1998) 325. 

5. C. B. Lang and P. Petreczky, Phys. Lett. 
387B (1996) 558. 

6. J. Jersak, C. B. Lang and T. Neuhaus, Phys. 
Rev. Lett. 77 (1996) 1933, Phys. Rev. D54 
(1996) 6909. 

7. J. Cox, T. Neuhaus and H. Pfeiffer, in prepa- 
ration 

8. A. M. Ferrenberg and R. H. Swendsen, Phys. 
Rev. Lett. 63 (1989) 1195. 

9. K. Binder, Rep. Prog. Phys. 50 (1987) 783. 

10. A. Billoire, T. Neuhaus and B. Berg, Nucl. 
Phys. B413 (1993) 795. 



